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(57) Abstract: A delayed decision feedback Log-MAP equalizer generates optimal soft information at the bit level for MPSK 
and MDPSK modulations over multipaih fading channels. The algorithms have a forward-backward recursion structure similar to 
conventional Log-MAP algorithms. To reduce complexity, the forward recursion incorporates the techniques of delayed decision 
feedback sequence estimation to construct a reduced state space. The backward recursion then takes this state space as given, in that 
it does not hypothesize feedbacks, and computes required soft information. 



BNSDOCID: <WO _0161952A1_I_> 



WO 01/61952 



PCT/US01/04312 



1 

DELAYED DECISION FEEDBACK LOG-MAP EQUALIZER 

BACKGROUND OF THE INVENTION 
The present invention relates generally to cellular and wireless communication 

and, more particularly, to an equalizer hypothesizing modulation symbols received over 
5 multipath fading channels- 

Signals in a wireless communication system are subject to a number of 
phenomena that degrade signal quality. Each signal is reflected from many different man-made 
and natural objects. The receiver thus receives a number of signals delayed by one or more 
signal periods, called "multipath", as each reflected signal is received. If the period of delay 

10 is more than the time required to transmit one symbol, producing intersymbol interference 

(ISI), then a receiver decoder may not be able to decode the symbol. This can cause poor 
sound or data quality to the user. Many different algorithms are used at the receiver to attempt 
to compensate for such effects. These techniques are discussed, for example, by G. L. Stuber 
in Principles of Mobile Communication, Chapter 6, Kluwer Academic Publishers, 1996. One 

1 5 such algorithm is the maximum likelihood sequence estimation (MLSE) or " Viterbi" algorithm. 

Digital cellular and personal communication systems based on IS- 1 36 and GSM 
(global system for mobile communications) require equalization to handle ISI arising from 
time dispersion. With the advent of the high date rates and the high level modulation in EDGE 
(enhanced data rates for global evolution) systems, the extent of ISI has increased 

20 considerably. To cover the same delay spread as today's GSM, EDGE receivers are specified 

to resolve a 5-tap channel. Furthermore, these systems are convolutional coded. To enable 
soft-decision decoding, it is, therefore, highly desirable for the equalizers to produce soft 
information at the bit level. 
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For M-ary modulation, a traditional MLSE equalizer requires A/- at each trellis 
stage in order to properly model an intersymbol interference (ISI) channel driven byL+\ M- 
ary phase-shift keying (MPSK) symbols. Because of the high computational and storage 
requirements this implies, the delayed decision feedback sequence estimation (DDFSE) 
equalizer uses a kind of hybrid of MLSE and DFE (decision feedback estimation) in which only 
the K + 1 most recent symbols are hypothesized as part of the state model, while the remaining 
L -K symbols required to compute each path metrics are determines from the surviving paths 
of each state. Therefore, a trellis of only A^' states is needed in this case and the L-K "DFE" 
symbols constitute side information of the states. DDFSE equalizers are discussed, for 
example, by A. Duel-Hallen and C Heegard in "Delayed Decision Feedback Sequence 
Estimation," IEEE Transactions on Communications, vol. 37, no. 5, pp. 428-436, May 1 989. 
One convenient way of describing the DDFSE operation is to say that the energy in the most 
recent K + 1 "MLSE" symbols is used to equalize the data, while the energy in the remaining 
"DFE" symbols is simply canceled. The performance of DDFSE equalizers can be improved 
by pre-filtering the received signals with minimum-phase filters to maximize the energy of the 
leading ISI taps. Such filtering techniques are discussed, for example, in copending 
Application Serial No. 09/378,314 entitled "Method and Apparatus for Computing Prefilter 
Coefficients" of K. C. Zangi et al. and filed Aug 20, 1999. 

In a convolutional coded system, aLog-MAP (maximum a posteriori) equalizer 
is preferred because of its ability to generate high-quality soft information for the Viterbi 
decoder. Traditional MAP algorithms, as described in L. R. Bahl, J. Cocke, F. Jelinek and J. 
Raviv (BCJR) in "Optimal Decoding of Linear Codes for Minimizing Symbol Error Rate," 
IEEE Transactions on Information Theory, vol. 20, pp. 284-287, Mar. 1974, are not very 
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suited to practical implementation because these algorithms require many multiplication and 
logarithm operations The Log-MAP algorithms circumvent this complexity problem by 
operating the signals in the logarithmic domain, as discussed, for example, by P. Robertson, 
E. Villebrun and P. Hoehner, in "A Comparison of Optimal and Sub-Optimal MAP decoding 
5 algorithms operating in the log domain," Proceedings of IEEE International Communications 

Conference '95, pp 1009-1013, June 1995. However, both approaches are designed to 
generate optimal soft information for trellis with static state spaces, i.e., the definitions of the 
states do not change over time. Moreover, these definitions are known a priori and are usually 
given exogenously. For example, the trellis of a convolutional code falls into this category: the 

10 states are defined by the contents of the encoding shift-registers and the complete trellis can 

be constructed even before the codewords are transmitted. For multipath equalization, the 
traditional MAP algorithms thus require the construction of the complete state space as in the 
case for MLSE equalizers. For many practical applications, this proves to be excessively 
complex. For example, a MAP equalizer for 8PSK with 5 channel taps (EDGE specifications) 

15 requires 4096 states. 

Since the Log-MAP equalizers suffer from the same complexity problems as 
the MLSE equalizers, it would be of great practical value if the principle of DDFSE could be 
combined with the Log-MAP algorithms to reduce complexity. An approach to this problem 
is discussed in the following. 

20 Several problems arise when this MLSE/DFE hybrid method is applied to a 

forward-backward Log-MAP equalizer. More specifically, the forward DDFSE recursion 
would use an anti-causal prefilter to create a minimum-phase channel response. ADDFSE run 
backwards through the trellis, on the other hand, would need a causal prefilter to crease a 
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maximum-phase channel, because in this case the MLSE symbols appear at the other end of 
the true state, as illustrated in Fig. ] . This results in two major issues: 

I The forward-backward Log-MAP algorithm breaks down. The two passes are 
now using two sets of data which are correlated with each other (because of 
5 the two prefilters with overlapping support). Because of this correlation, the 

Markovian property of the data sequence, which means that the future data 
symbols are independent of the past data symbols given the current state, no 
longer holds. This Markovian property plays an essential role in the derivation 
of the forward-backward algorithm, without which the Log-MAP algorithm is 
10 invalid. 

2. The state space from the two recursions do not match. A complete state space 
is used in a Log-MAP equalizer, which guarantees the forward and backward 
recursions work on the same state space. On the other hand, for a naive 
delayed decision feedback Log-MAP equalizer, only the "MLSE" symbols can 
1 5 be matched up. The "DFE" symbols setup by the forward and backward 

recursions generally do not agree with each other. The overlapping supports 
of the two prefilters makes the situation even more confusing. 

The definitions of the states in a DDFSE trellis are constructed endogenously 
by the equalization recursion itself and, therefore, become unknown a priori. The traditional 
20 MAP algorithms are thus not directly applicable in these situations. In order to extract soft 

information, ad hoc post-processing methods such as the soft-output Viterbi algorithms 
(SOVA) are usually used instead The SOVA algorithm is discussed, for example, by J 
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Jagenauer and P. Hoehner in "A Viterbi Algorithm with Soft-Decision Outputs and its 
Applications/' Proceedings of IEEE Global Telecommunications Conference '89, pp. 1680- 
1 686, Nov. 1 989. These algorithms, however, are sub-optimal and their shortcomings become 
more apparent in iterative processing settings. 
5 Accordingly, there is a need for an equalizer and method that generates optimal 

soft information at the bit level for MPSK and MDPSK modulations over multipath fading 
channels. 

SUMMARY OF THE INVENTION 
The present invention meets this need by providing a delayed decision feedback 
10 Log-MAP equalizer that generates optimal soft information at the bit level for MPSK and 

MDPSK modulations over multipath fading channels. 

Broadly, there is disclosed in accordance with one aspect of the invention a 
method of hypothesizing modulation symbols received over multipath fading channels. An 
estimate of a received sequence of symbols is calculated and stored using stored channel 
1 5 coefficients and a previous state path history of modulation symbols. For each possible current 

state symbol, a branch metric is calculated and stored from a prior state to a current state using 
the calculated estimate of the received sequence of symbols. A forward recursion is performed 
by calculating a forward state metric using the branch metric from the prior state to the current 
state. A backward recursion is performed by calculating a backward state metric using the 
20 branch metric from the prior state to the current state. 

It is a feature of the invention that the estimate of a received sequence of 
symbols calculated and stored using stored channel coefficients and a previous state path 
history of modulation symbols is implemented using a delayed decision feedback sequence 
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estimation equalizer. 

It is another feature of the invention that the estimate of a received sequence 
of symbols calculated and stored using stored channel coefficients and a previous state path 
history of modulation symbols is implemented using a maximum likelihood sequence estimation 
equalizer for most recent symbols in the sequence of symbols and using a decision feedback 
equalizer for remaining symbols in the sequence of symbols. 

It is a further feature of the invention that a Log-MAP equalizer is used to 
calculate and store the branch metric, the forward state metric and the backward state metric. 

It is yet another feature of the invention that a posteriori log-likelihood ratio 
of each bit in the hypothesized modulation symbols is calculated using the calculated forward 
state metric and the calculated backward state metric. 

It is still another feature of the invention that the forward recursion is executed 
prior to the backward recursion. 

It is still an additional feature of the invention that the backward recursion is 
executed prior to the forward recursion. 

It is yet still another feature of the invention that the previous state path history 
of modulation symbols is updated using the calculated branch metric from a prior state to a 
current state using the calculated estimate of the received sequence of symbols. 

It is still a further feature of the invention that the modulation symbols received 
over multipath fading channels comprise M-ary modulated symbols. 

It is yet a further feature of the invention that the modulation symbols received 
over multipath fading channels comprise M-ary differential PSK modulated symbols. 

There is disclosed in accordance with another aspect of the invention, in a 
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receiver receiving modulation symbols transmitted over multipath fading channels, an equalizer 
including a programmed controller. The programmed controller calculates and stores an 
estimate of a received sequence of symbols using stored channel coefficients and a previous 
state path history of modulation symbols A branch metric is calculated and stored for each 
5 possible current state symbol from a prior state to a current state using the calculated estimate 

of the received sequence of symbols. A forward recursion is performed by calculating a 
forward state metric using the branch metric from the prior state to the current state. A 
backward recursion is performed by calculating a backward state metric using the branch 
metric from the prior state to the current state. 

1 0 There is disclosed in accordance with a further aspect of the invention a method 

of generating soft information at bit level of modulation symbols received over multipath 
fading channels. An estimate of a received sequence of symbols is calculated and stored using 
stored channel coefficients and a previous state path history of modulation symbols, each 
symbol comprising an M-ary modulated symbol representing select bits c. A branch metric for 

1 5 each possible current state symbol is calculated and stored from a prior state to a current state 

using the calculated estimate of the received sequence of symbols. A forward recursion is 
performed by calculating a forward state metric using the branch metric from the prior state 
to the current state. A backward recursion is performed by calculating a backward state metric 
using the branch metric from the prior state to the current state. A posteriori log-likelihood 

20 ratio of each select bit c in the estimated modulation symbols is calculated using the calculated 

forward state metric and the calculated backward state metric. 

Particularly, a DDFSE method is incorporated with a forward recursion to 
construct a reduced state space Starting with a complete state space, the M 1 states are 
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partitioned into h/F subsets. Each of the subsets consists of the M L " K) states with the same K 
leading hypothesized symbols while the remaining L-K hypotheses are completely different. 
Since the minimum-phase prefilter shifts most of the energy to the leading taps, the last L-K 
hypotheses have little effect on the recursion metrics. To reduce computation, L -K decision 
feedbacks are then used to select one representative state for each subset. 

The same reduced state space is used by the backward recursion. Using this 
structure, the decision feedback side information for the backward recursion is not associated 
with the surviving paths in the backward trellis, but rather is determined by the forward 
recursion to represent backward prediction symbols along each path. That is, as shown in Fig. 
2, the DFE symbols appear on the left of the MLSE symbols for both forward and backward 
recursions, meaning that a minimum-phase channel is appropriate for both directions. 
Therefore, both directions can use the same prefilter and the two concerns described in the 
previous section no longer apply. 

The disclosed invention generalizes the traditional BCJR paradigm to this 
reduced-state trellis, which is referred to as a DDF-Log-MAP (delayed decision feedback Log- 
MAP) algorithm. The disclosed method has a forward-backward recursion structure similar 
to the classic Log-MAP algorithms. To reduce complexity, however, the forward recursion 
incorporates DDFSE to construct a reduced state space. The backward recursion then takes 
this state space as given, i.e., it does not hypothesize feedbacks, and computes the required soft 
information. The description includes exemplary applications of the invention to MPSK and 
M-ary differential phase-shift keying (MDPSK) modulations, respectively, with a reduced state 
space equivalent to that of a DDFSE with 2 MLSE taps. Generalizations to other modulations 
and other choices of reduced state spaces are straightforward. 
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Further features and advantages of the invention will be readily apparent from 
the specification and the drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates state configurations for forward and backward recursions using 
a naive approach to combining DDFSE and Log-MAP algorithms; 

Fig. 2 illustrates state configurations for forward and backward recursions using 
an approach to combining DDFSE and Log-MAP algorithms in accordance with the invention; 

Fig. 3 is a block diagram of a mobile terminal implementing the method 
according to the invention; 

Fig. 4 is a block diagram of a system model for a wireless communication 
system implementing the method according to the invention; 

Fig. 5 illustrates trellis notation used in explaining the method according to the 

invention; and 

Fig. 6 is a flow diagram illustration of a program run by the processor of Fig. 
3 to implement the method according to the invention. 

DETAILED DESCRIPTION OF THE INVENTION 
In accordance with the invention, a DDFSE method is incorporated with the 
forward recursion to construct a reduced state space. Starting with a complete state space, 
theA^ states are partitioned into subsets. Each of the subset consists of the states 
with the same K leading hypothesized symbols while the remaining L-K hypotheses would be 
completely different. Since the minimum-phase prefilter shifts most of the energy to the 
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leading taps, the last L - K hypothesis have little effect on the recursion metrics. To reduce 
computation, L-K decision feedbacks are then used to select one representative state for each 
subset. 

The same reduced state space is used by the backward recursion. The decision 
5 feedback side information for the backward recursion is not associated with the surviving paths 

in the backward trellis, but rather is determined by the forward recursion to represent 
backward prediction symbols along each path. That is, as shown in Fig. 2, the DFE symbols 
appear on the left of the MLSE symbols for both forward and backward recursions, meaning 
that a minimum-phase channel is appropriate for both directions. Therefore, both directions 

10 can use the same prefilter. 

Fig. 3 is a block diagram of a typical mobile terminal 1 0. The mobile terminal 
10 includes an antenna 12, a receiver 16, a transmitter, 18, a speaker 20, a processor 22, a 
memory 24, a user interface 26 and a microphone 32. The antenna is configured to send and 
receive radio signals between the mobile terminal 1 0 and a wireless network (not shown). The 

1 5 antenna 1 2 is connected to a duplex filter 1 4 which enables the receiver 1 6 and the transmitter 

18 to receive and broadcast (respectively) on the same antenna 12. The receiver 16, in 
conjunction with the processor 22, demodulates, demultiplexes and decodes the radio signals 
into one or more channels. Such channels include a control channel and a traffic channel for 
speech or data. The speech or data are delivered to the speaker 20 (or other output device, 

20 such as a modem or fax connector). 

The receiver 16 delivers messages from the control channel to the processor 
22. The processor 22 controls and coordinates the functioning of the mobile terminal 10 and 
is responsive to messages on the control channel using programs and data stored in the 
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memory 24, so that the mobile terminal 10 can operate within the wireless network (not 
shown). The processor 22 also controls the operation of the mobile terminal 10 and is 
responsive to input from the user interface 26. The user interface 26 includes a keypad 28 as 
a user-input device and a display 30 to give the user information. Other devices are frequently 
5 included in the user interface 26, such as lights and special purpose buttons The processor 22 

controls the operations of the transmitter 1 8 and the receiver 1 6 over control lines 34 and 36, 
respectively, responsive to control messages and user input. 

The microphone 32 (or other data input device) receives speech signal input and 
converts the input into analog electrical signals. The analog electrical signals are delivered to 

10 the transmitter 18. The transmitter 18 converts the analog electrical signals into digital data, 

encodes the data with error detection and correction information and multiplexes this data with 
control messages from the processor 22. The transmitter 18 modulates this combined data 
stream and broadcasts the resultant radio signals to the wireless network through the duplex 
filter 14 and the antenna 12. 

15 While the present invention is described relative to a cellular communication 

system, the invention is applicable to other wireless communication systems having 
communication channels subject to intersymbol interference (ISI). 

Digital cellular and personal communication systems based on IS- 1 3 6 and GSM 
(global system for mobile communications) require equalization to handle ISI arising from 

20 time dispersion. With the advent of the high date rates and the high level modulation in EDGE 

(enhanced data rates for global evolution) systems, the extent of ISI has increased 
considerably. To cover the same delay spread as today's GSIvl, EDGE receivers are specified 
to resolve a 5-tap channel. Furthermore, these systems are convolutional coded, as illustrated 



■■NSDOCID: <W<D 



_0161952A1J_> 



WO 01/61952 PCT/US01/04312 

12 

in Fig. 4. To enable soft-decision decoding, it is, therefore, highly desirable for the equalizers 
to produce soft information at the bit level. 

As illustrated in the system mode! of Fig. 4, data is input to a convolutional 
encoder 40. The encoded data is output through a channel interleaver 42 to a modulator 44. 
The modulator 44 modulates the bits in a conventional manner to produce a signal x n to be 
transmitted via a multipath fading channel represented at 46. A received sequence /*„ is input 
to an equalizer 48 in accordance with the invention to produce the soft information. The 
output of the equalizer 48 is delivered through a channel de-interleaver 50 to a convolutional 
decoder 52. 

For // = 1, 2, . . N, the received sequence r n is related to the transmitted 

sequence x n by 

L 
/=0 



where ^/=0J are complex-valued multipath fading coefficients and w n is complex- 
valued, zero-mean white Gaussian noise. As discussed above, to improve performance, 
15 minimum-phase pre-filtering can be used to maximize the energy in the leading coefficients, 
i. e., the energy in h l for small /. 

For MPSK modulation, where M = 2\ the transmitted signal x n is uniquely 

2n 

specified by the input (/-tuple, (c n]y c n2 , .,c nef ) e {0,1}", via x n = m s = e~^ S , where 

S = q *2 <M 4- C 2 *2*~ 2 -f ...+ c q . A reverse <y-ary to binary mapping is defined as£ 6 (s) = 
20 c h ,/>= 1,2 q. 
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For MDPSK modulation, where M = 2 9 , the transmitted signal x n is given by 

X n = U nJ (2) 

where v n is an MPSK symbol uniquely specified by the input <?-tuple, (c„ l5 c n2 , c^) e 

{0,1 } q , via X n = m s = e A/ * , where S = c x *2 g ~ ] + c 2 *2 9 " 2 + ...+ £* 9 . A reverse ^-ary 

to binary mapping is defined as B b (s) = c b , b = 1, 2, . . q. 

The problem is to find, for all // and b, the a posteriori log-likelihood ratio of 

the bit c nh . 

A(c n<> )Alog-J -T-4- 

given the complete received sequence = (r,, r 2 , . . r N ) and the a priori log-likelihood 



(3) 



1 0 ratios 



Using the traditional log-MAP algorithm, referring to the notations described 
in Fig. 4, the forward state metric (5), backward state metric (6), and the branch metric (7, 
8) are defined as 

15 a „(j) = \ogPr {S n = j\ r "} (5) 

A^>iogPr{rlk = y} («) 
r n (ij)= io g Pr{5 n = y,A-„|s„., = /} (7) 

= log Pr {S„ = j 1 = /} + log Pr {r„ | S„_, = /, 5„ = 7} , (8) 
where f" m represents the sequence (/„„ /„, + „ . . ., / „). The forward and backward recursions 

20 are: 

«„0)= COM(a „_,(/)+ r „VJ)) M»= 1,2,... ,N (9) 
/?„_,(/) = COM{fi.U)+ r„VJ)), Mr>= N,N-1,...,\ (10) 
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where the soft combining operation COM is defined as 

COA/(*,)Alog(5>*') . (11) 
Efficient implementation of this operation is given using: 

COM(x 9 y) = rnax(.x 5 j/)+ log(l + e~ lx ~ yl ) . (12) 

5 That is, the soft combiner first finds the maximum of the two inputs, as in the case for the 

Viterbi algorithm, and then adds a correction term. Many efficient implementations of this 
operation are discussed, for example, by P. Robertson, E. Villebrun and P. Hoehner, in "A 
Comparison of Optimal and Sub-Optimal MAP Decoding Algorithms Operating in the Log 
Domain," Proceedings of IEEE International Communications Conference l 95, pp. 1009- 
10 1013, June 1995. One ofthe approach ignores the correction term completely to obtain a very 

computationally efficient algorithm, which is widely referred to as the Log-Max algorithm. 
The invention disclosed herein can be used in conjunction with any of these techniques. 

If the trellis begins at known state B and ends at known state E, the two 
recursions are initialized by 

foe, for j '= B foe , for j = E 

oWy [0, otherwise HnKJJ [0, otherwise v ' 

where °° represents a very large number permitted by the processor 22 or a special symbol 
recognizable by the processor 22 The a posteriori state log-likelihood and the a posteriori 
edge log-likelihood from state / to state j at time n can be then computed as 

iogPr{s„ = yjrf} = ff.cm AC/) (") 
io g Pr{5„_, = i,s n = j\ r ") = «„_,(/) + r.(/,y) + AO") • 0 5 ) 

The following describes the use of DDF-Log-M AP for MPSK in accordance 



20 
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with the invention. 

Under the structure of DDFSE with 2 MLSE taps, there are M states at each 
stage of the DDFSE trellis. Associated with each state / = 0, 1, . . . , M - 1, at the (?i - l)-th 
stage is a path history of L-tuple MPSK symbols corresponding to hypothesized transmitted 
5 signals from stage // - L to // - 1 : 

P t (n- 1)= [x,(n- L), Xi (n- L+ 1). ..,*,(«- 1) = m] . (16) 

For branch metrics, since all the paths that lead to the same end state j are driven by the same 
hypothesized symbol m r the first term of equation (8) is determined by: 

logPr[s„ = /IS,,., = /}= logPr{*„ = m,} 

= £logPr{ c - 

= I B b (J)\og U I logPr {c nb = 0} 

10 where V(c nA ) is the a priori log-likelihood ratio for the bit c nb and B b (J) is the #-ary to binary 

function defined above. For this equalizer, the second term of equation (8) can be written as: 



L 

logPr {r n \S n . { = i,S„ = j} = K' - r n - h,m } - - I) 

= Ar-|r„(/)| 2 + 2Re{/z>;/;(/)} , 



(18) 
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where Re means taking the real part and fjij is defined as 

= r n -Y< V/("-0 (19) 

which is the estimated received signal after decision-feedback cancellation. The two terms, 
2/= A*/ ( w ~ 0 an d 2R e {^o^J^iO)} > can be efficiently computed by, for example, 

5 techniques discussed in copending Application Serial No. , entitled "Methods, 

Receivers and Equalizers Having Increased Computation Efficiency" to K C. Zangi et al. and 
filed Feb. 8, 2000 (Attorney Docket No. 8 1 94-354), which is incorporated by reference herein 
For the forward recursion, since COM{a + jc.) = a + COM(x i ) , the«-th 

stage of the forward recursion can then be re-arranged as 

10 a,XJ)=COM(a„_ l (i)- \r„(if + 2 Re{/;»; (/)})+ £ AC/)K(c*). (20) 

For the decision feedback update, the new path history for state j is defined as 
P J {n)=[x s (n- L+ l),...,jr f (»- 2),ni„m / ] , 

(2D 

where s is the index of the prior state from which the surviving branch is originated and: 

15 s= ,=o m . a 5-,"' (*~.<o- \ p M 2 + 2Re{«»;^(o}) • 

(22) 

That is, the surviving branch into state/ is defined as the branch ending at J with the maximum 
accumulated metric 

(*„-,(')- fc(0f + 2 Re{/7»;(/)}) . (23) 
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Since finding the maximum is part of the COM operation, the last expression does not incur 
any new complexity. 

For the backward recursion, the (/?-l)-th stage can be re-arranged as: 

A-,(0 = CoUfi.Ul+t. **C/)V(0+ 2Re{/*>X(0}j - |r„(/)| 2 . (24) 

5 That is, the backward recursion does not hypothesize feedbacks and instead re-uses the 

modified received sequence P n (i) computed by the forward recursion. 

For the soft output, the a posteriori log-likelihood ratio of the bit c nb is then 

given by 

A(c„,)= COM (*„(*)+ /?„(*))- COM (a (25) 

10 That is, for each 6=1, 2,...,tf, the a posteriori state log-likelihoods (a n (s) + fi n (5)), , for 

s = 0, 1 , . . M - 1 , are first separated into two sets such that one set has the Z>-th position in the 
binary representation of 5 being 1 and the other set has that position being 0. These two sets 
of log-likelihoods are then processed by the soft combiner separtely. Finally, the a posteriori 
log-likelihood for bit Cnb is obtained by taking the difference of the two combined results. The 
15 extrinsic information for the bit c nb can be obtained as 

W{c nh )= A(0- V(0. (26) 

The above described method is formulated to execute the forward recursion 
first and then backward recursion in the above. It is straightforward to generalize it to execute 
the backward recursion first and followed by the forward recursion. 
20 Using proper implementations for the COM operations and techniques to 

compute 
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htXijn- /) and 2Re^*/;?*r w (/) , the disclosd method is comparable to the 

traditional or SOVA-type DDFSE in total computational complexity. However, in addition 
to those variables generally used in MLSE or DDFSE equalizers, the disclosed method stores 
the forward state metrics a n (JJ if the forward recursion is executed before the backward 

5 recursion. If the backward recursion is executed before the forward recursion, then the 

backward state metrics (5 n (i) are stored. For systems with small blocks, e. g., D-AMPS, this 
storage requirement is not too excessive. For systems with large blocks, the disclosed method 
can be applied in a sliding-window style as discussed by, for example, S. Benedetto, G. 
Montorsi, D Divsalar and F. Pollara in "Soft-Output Decoding Algorithms in Iterative 

1 0 Decoding of Turbo Codes," The Telecommunications and Data A cquisitioti Progress Report, 

Jet Propulsion Laboratory, California Institute of Technology, vol. 42-124, pp. 63-87, Feb. 
1996, to reduce the required strorage size. 

The following describes the use of DDF-Log-MAP for MDPSK in accordance 
with the invention 

1 5 Under the structure of DDFSE with 2 MLSE taps, there are M states at each 

stage of the DDFSE trellis. Associated with each state / = 0, 1, . . , M - 1, at the (n - l)-th 
stage is a path history of L-tupie MPSK symbols corresponding to hypothesized transmitted 
signals from stage n - L to // - 1 

P t (n- 1) = [jc,(>7- L\x,(n- L+ 1).. 2),x,(w- 1) = m] . (27) 
20 For branch metrics, because of the differential encoding, the transition from 

state / to state / is driven by //„ = m 6 , where 8 € {0,1, . . . , M - 1 } is obtained by 

8 = ('-./XwA('- 7)modM. (28) 
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where mod is the modulo operation. Therefore, the first term of equation (8) can be computed 
as 

logPr{s„ = y|s„., = /}= logPr{w„ = m s ) 

= ZlogPr{c„, = ^((/-y)^)} 

b=) V 3 

h=\ 

where V(c nb ) is the a priori log-likelihood ratio for the bit and B b is the <?-ary to binary 

5 function defined above. Note that ^ A B h i(i — j) )K(c n6 )can only assume one of the 

b=\ 

following M values: 

2>„(<?)V(c„„)= X v (0, * = 0,1,...,M-1 (30) 



which can be pre-computed using only 



M 

((7-2) — + ] 



real additions. 



10 



Since all paths that lead to the same end state j hypothesize the same 
transmitted symbol m p for this equalizer, the second term of equation (8) can be computed as: 



IogPi^^,,., = i,S„ = j} = K> - r„ - - £ V,(»" /) 

= ^-|/;,(/)| 2 + 2Re{A»;(/)}, 
where Re means taking the real part and V is defined as 

L 

r„(')=/;-Zv,(»-/) 

/=! 



(31) 



(32) 
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which is the estimated received signal after decision-feedback cancellation. The two terms, 
~ O and 2Re|/?o/i7jr fI (o} . can be efficiently computed by, for example, 

techniques discussed in the aforementioned reference incorporated herein. 

For the forward recursion, the /Mh stage can then be formulated as 

a.(J)' COM[ (*„_,(/)- |r.(0f) + 2Re{/ J >7„(,)) + £ B b ((i- j) M )v(c nb ^3) 



For the decision feedback update, the new path history for state; is defined as 
/>(/?)= [x,(/7- L + 1) 2),m J> m.] , (34) 

where 5 is the index of the prior state from which the surviving branch is originated and: 

S= , = o m .^-r{("»-' ( 'H^ • 05) 

1 0 That is, the surviving branch into state j is defined as the branch ending at j with the maximum 

accumulated metric 



((".-.CO - k(0| 2 ) + 2Re{/7>X(/)} + X 2^«/- y), y )F(c w6 )) . 



(36) 



Since finding the maximum is part of the COM operation, the last expression does not incur 
1 5 any new complexity. 

For the backward recursion, since COM {a + x. ) = <7 + COM(x { ) ,the(»- 
l)-th stage of the backward recursion can be rearranged as: 
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£,-,(')= COM^ (37) 

That is, the backward recursion does not hypothesize feedbacks and instead re-uses the 
modified received sequence P n (i) computed by the forward recursion. 

For the soft output, the a posteriori log-likelihood ratio of the bit c nh is then 



5 given by 



A(0= COM (a n _ } (0+rAiJ)+fi„U)) 



com {a^o) + rAU)+fi.U)) 



(38) 



where yJJJ) as defined in equation (8) is simply the sum of equations (29) and (31) That is, 
the A/ 3 a posteriori edge log-likelihoods can be computed as («„_,(/)+ y n (i\j) + PJJ)), for each 
l : = 0, 1, A/ - 1 and / = 0, 1, A/ - 1. For each 6 = 1, 2, these log-likelihoods are first 
10 separated into two sets such that one set has the Z>-th position in the binary representation 

of (j — y ),U being 1 and the other set has that position being 0. These two sets of log- 
likelihoods are then processed by the soft combiner separately. Finally, the a posteriori log- 
likelihood for bit Cnb is obtained by taking the difference of the two combined results. 

It should be noted that the computational complexity of equation (38) can 
15 actually be reduced by first computing the a posteriori log-likelihoods of the driving symbol 

u n and then deriving the bit log-likelihoods from these symbol log-likelihoods. For M= 8, this 
two-step approach can reduce the complexity of equation (38) by 60%. 

This method is demonstrated in the following. First, the M driving symbol log- 
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likelihoods are obtained as: 

r„(*)= c5Mk,('>r»(',('>4)+ *>„)), m 



for 6 = 0, ],..., A/-1 . The a posteriori log-likelihood for bit c nb can then be computed using 

a (O = COM(r„(^))- C<9M(r„ (<?))• (40) 

5 That is, for each 6=1,2,.., <:/, these symbol log-likelihoods are first separated into two sets 

such that one set has the A-th position in the binary representation of 6 being 1 and the other 
set has that position being 0. These two sets of log-likelihoods are then processed by the soft 
combiner separately. Finally, the a posteriori log-likelihood for bit c nb is obtained by taking 
the difference of the two combined results. 

10 The extrinsic value for bit c nb can be obtained as 

Similar to the DDF-Log-M AP method for MPSK equalization discussed above, 
the disclosed method for MDPSK equalization can also be generalized to execute the 
backward recursion first and followed by the forward recursion. Similar comments on 
15 complexity from that discussion apply to the MDPSK equalization method as well. 

Referring to Fig. 6, a flow diagram illustrates a program implemented in the 
processor 22 of Fig. 3 for implementing the method of delayed decision feedback Log-MAP 
equalization in accordance with the invention. As will be apparent, the equalization method 
is done in coordination with hardware in the receiver 16 of Fig. 3. In describing the method, 
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reference will be made to equations for both MPSK and, in parenthesis, MDPSK 

The program begins at a start node 60 when a signal is received. The program 
has a first loop 62 performing the forward recursion for n=l, 2, ...N and, for each n, j=0, 1, 
...M-l. A second loop 64 performs the backward recursion for n=N, N-l, ...1 and, for each 
5 n> i=0, 1, ...M-l. 

The first loop 62 begins at a block 66 that calculates and stores the received 
signal estimate using equation 1 9 (32) based on the path history from equation 16 (27). Then, 
at a block 68, the branch metric is calculated using equations 1 7 and 1 8 (29 and 3 1 ). Finally, 
at a block 70, the forward recursion is performed. The forward recursion implements equation 
10 20 (33) and then for each state j updates the path history of equation 21 (34) using equation 

22 (35). A decision block 72 loops back to the block 66 until the loop 62 completes a forward 
trellis for all n and j and the estimated received signal values, forward metrics and path histories 
are stored in memory 24, see Fig. 3, for each state and time. Thereafter, the second loop 64 
is performed to run a backward trellis. 
1 5 The second loop 64 begins at a block 74 that performs the backward recursion. 

The backward recursion implements equation 24 (37) for each state i. Then, at a block 76, the 
soft output is calculated by determining the a posteriori log likelihood ratio using equation 25 
(38 or 39 and 40). Finally, extrinsic information is obtained for each bit at a block 78 using 
equation 26 (41). A decision block 80 loops back to the block 74 until the loop 64 completes 
20 the backward trellis for all n and i. The routine is then done. 

As is apparent, to implement the backward recursion first the functions of 
blocks 70 and 74 would be reversed. 

Thus, in accordance with the invention, an efficient method generates optimal 
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soft information at the bit level for MPSK and MDPSK modulations over multipath fading 
channels. The algorithms have a forward-backward recursion structure similar to conventional 
Log-MAP algorithms To reduce complexity, the forward recursion incorporates the 
techniques of DDFSE to construct a reduced state space. The backward recursion then takes 
this state space as given, in that it does not hypothesize feedbacks, and computes required soft 
information. 
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CLAIMS 

1 . A method ofhypothesizing modulation symbols received over multipath 
fading channels, comprising: 

calculating and storing an estimate of a received sequence of symbols using 
stored channel coefficients and a previous state path history of modulation symbols; 

for each possible current state symbol, calculating and storing a branch metric 
from a prior state to a current state using the calculated estimate of the received sequence of 
symbols, 

performing a forward recursion by calculating a forward state metric using the 
branch metric from the prior state to the current state; and 

performing a backward recursion by calculating a backward state metric using 
the branch metric from the prior state to the current state. 

2. The method of claim 1 wherein calculating and storing the estimate of 
a received sequence of symbols using stored channel coefficients and a previous state path 
history of modulation symbols is implemented using a delayed decision feedback sequence 
estimation equalizer. 

3. The method of claim 1 wherein calculating and storing the estimate of 
a received sequence of symbols using stored channel coefficients and a previous state path 
history of modulation symbols is implemented using a maximum likelihood sequence estimation 
equalizer for most recent symbols in the sequence of symbols and using a decision feedback 
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equalizer for remaining symbols in the sequence of symbols. 

4 . The method of claim 1 wherein calculating and storing the branch metric 
from a prior state to a current state using the calculated estimate of the received sequence of 
symbols, calculating the forward state metric using the branch metric from the prior state to 
the current state, and calculating the backward state metric using the branch metric from the 
prior state to the current state are implemented using a Log-MAP equalizer. 

5. The method of claim ) further comprising calculating a posteriori log- 
likelihood ratio of each bit in the estimated modulation symbols using the calculated forward 
state metric and the calculated backward state metric. 

6. The method of claim 1 wherein performing the forward recursion by 
calculating the forward state metric using the branch metric from the prior state to the current 
state is executed prior to performing the backward recursion by calculating the backward state 
metric using the branch metric from the prior state to the current state. 

7. The method of claim 1 wherein performing the backward recursion by 
calculating the backward state metric using the branch metric from the prior state to the 
current state is executed prior to performing the forward recursion by calculating the forward 
state metric using the branch metric from the prior state to the current state. 

8. The method of claim 1 further comprising updating the previous state 
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path history of modulation symbols using the calculated branch metric from a prior state to a 
current state using the calculated estimate of the received sequence of symbols. 

9. The method of claim 1 wherein the modulation symbols received over 
multipath fading channels comprise M-ary modulated symbols. 

5 10. The method of claim 1 wherein the modulation symbols received over 

multipath fading channels comprise M-ary differential PSK modulated symbols 

11. In a receiver receiving modulation symbols transmitted over multipath 
fading channels, an equalizer comprising: 

a controller programmed to: calculate and store an estimate of a received 
10 sequence of symbols using stored channel coefficients and a previous state path history of 

modulation symbols; for each possible current state symbol calculate and store a branch metric 
from a prior state to a current state using the calculated estimate of the received sequence of 
symbols; perform a forward recursion by calculating a forward state metric using the branch 
metric from the prior state to the current state; and perform a backward recursion by 
1 5 calculating a backward state metric using the branch metric from the prior state to the current 

state. 

12. The equalizer of claim 1 1 wherein the controller calculates the estimate 
of a received sequence of symbols using stored channel coefficients and a previous state path 
history of modulation symbols using a delayed decision feedback sequence estimation. 
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1 3 . The equalizer of claim 1 1 wherein the controller calculates the estimate 
of a received sequence of symbols using stored channel coefficients and a previous state path 
history of modulation symbols using maximum likelihood sequence estimation for most recent 
symbols in the sequence of symbols and using a decision feedback equalization for remaining 
symbols in the sequence of symbols. 

14. The equalizer of claim 1 1 wherein the controller calculates the branch 
metric from a prior state to a current state using the calculated estimate of the received 
sequence of symbols, calculates the forward state metric using the branch metric from the prior 
state to the current state, and calculates the backward state metric using the branch metric 
from the prior state to the current state using Log-MAP equalization. 

1 5 The equalizer of claim 1 1 wherein the controller calculates a posteriori 
log-likelihood ratio of each bit in the estimated modulation symbols using the calculated 
forward state metric and the calculated backward state metric. 

1 6. The equalizer of claim J 1 wherein the controller performs the forward 
recursion by calculating the forward state metric using the branch metric from the prior state 
to the current state prior to performing the backward recursion by calculating the backward 
state metric using the branch metric from the prior state to the current state. 
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1 7. The equalizer of claim 1 1 wherein the controller performs the backward 
recursion by calculating the backward state metric using the branch metric from the prior state 
to the current state prior to performing the forward recursion by calculating the forward state 
metric using the branch metric from the prior state to the current state. 



5 1 8. The equalizer of claim 1 1 wherein the controller updates the previous 

state path history of modulation symbols using the calculated branch metric from a prior state 
to a current state using the calculated estimate of the received sequence of symbols 



19. The equalizer of claim 1 1 wherein the modulation symbols received over 
multipath fading channels comprise M-ary modulated symbols. 

10 20. The equalizer of claim 1 1 wherein the modulation symbols received over 

multipath fading channels comprise M-ary differential PSK modulated symbols. 



21. A method of generating soft information at bit level of modulation 
symbols received over multipath fading channels, comprising: 

calculating and storing an estimate of a received sequence of symbols using 
15 stored channel coefficients and a previous state path history of modulation symbols, each 

symbol comprising an M-ary modulated symbol representing select bits c; 

for each possible current state symbol, calculating and storing a branch metric 
from a prior state to a current state using the calculated estimate of the received sequence of 
symbols; 



3NSDOCID: <WO 0161952A1_I_> 



WO 01/61952 PCT/US01/04312 

30 

performing a forward recursion by calculating a forward state metric using the 
branch metric from the prior state to the current state; 

performing a backward recursion by calculating a backward state metric using 
the branch metric from the prior state to the current state; and 
5 calculating a posteriori log-likelihood ratio of each select bit c in the estimated 

received modulation symbols using the calculated forward state metric and the calculated 
backward state metric. 

22. The method of claim 2 1 wherein calculating and storing the estimate of 
a received sequence of symbols using stored channel coefficients and a previous state path 

10 history of modulation symbols is implemented using a delayed decision feedback sequence 

estimation equalizer. 

23 . The method of claim 2 1 wherein calculating and storing the estimate of 
a received sequence of symbols using stored channel coefficients and a previous state path 
history of modulation symbols is implemented using a maximum likelihood sequence estimation 

15 equalizer for most recent symbols in the sequence of symbols and using a decision feedback 

equalizer for remaining symbols in the sequence of symbols. 

24. The method of claim 21 wherein calculating and storing the branch 
metric from a prior state to a current state using the calculated estimate of the received 
sequence of symbols, calculating the forward state metric using the branch metric from the 

20 prior state to the current state, and calculating the backward state metric using the branch 



BNSDOCID: <WO 0161952A1 J_> 



WO 01/61952 PCTYUS01/04312 

31 

metric from the prior state to the current state are implemented using a Log-MAP equalizer 

25. The method of claim 21 wherein performing the forward recursion by 
calculating the forward state metric using the branch metric from the prior state to the current 
state is executed prior to performing the backward recursion by calculating the backward state 
metric using the branch metric from the prior state to the current state. 

26. The method of claim 2 1 wherein performing the backward recursion by 
calculating the backward state metric using the branch metric from the prior state to the 
current state is executed prior to performing the forward recursion by calculating the forward 
state metric using the branch metric from the prior state to the current state. 

27. The method of claim 2 1 further comprising updating the previous state 
path history of modulat ion symbols using the calculated branch metric from a prior state to a 
current state using the calculated estimate of the received sequence of symbols. 

28 . The method of claim 2 1 wherein the modulation symbols received over 
multipath fading channels comprise M-ary differential PSK modulated symbols. 
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29. A mobile terminal receiving modulation symbols transmitted over 
rnultipath fading channels, comprising: 

a receiver for receiving modulation symbols transmitted over rnultipath fading 
channels and including an equalizer calculating and storing an estimate of a received sequence 
of symbols using stored channel coefficients and a previous state path history of modulation 
symbols, calculating and storing a branch metric from a prior state to a current state using the 
calculated estimate of the received sequence of symbols for each possible current state symbol; 
performing a forward recursion by calculating a forward state metric using the branch metric 
from the prior state to the current state; and performing a backward recursion by calculating 
a backward state metric using the branch metric from the prior state to the current state. 

30. The mobile terminal of claim 29 wherein the equalizer calculates the 
estimate of a received sequence of symbols using stored channel coefficients and a previous 
state path history of modulation symbols using a delayed decision feedback sequence 
estimation 



3 1 The mobile terminal of claim 29 wherein the equalizer calculates the 
estimate of a received sequence of symbols using stored channel coefficients and a previous 
state path history of modulation symbols using maximum likelihood sequence estimation for 
most recent symbols in the sequence of symbols and using a decision feedback equalization for 
remaining symbols in the sequence of symbols. 

32. The mobile terminal of claim 29 wherein the equalizer calculates the 
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branch metric from a prior state to a current state using the calculated estimate of the received 
sequence of symbols, calculates the forward state metric using the branch metric from the prior 
state to the current state, and calculates the backward state metric using the branch metric 
from the prior state to the current state using Log-MAP equalization 

5 33. The mobile terminal of claim 29 wherein the controller calculates a 

posteriori log-likelihood ratio of each bit in the hypothesized modulation symbols using the 
calculated forward state metric and the calculated backward state metric. 

34. A method of generating soft information at bit level of modulation 

symbols received over multipath fading channels, comprising; 
10 performing a forward recursion by calculating and storing an estimate of 

received symbols using decision feedback cancellation, each symbol comprising an M-ary 

modulated symbol representing select bits c, calculating and storing branch metrics using the 

calculated estimate of the received symbols, and calculating forward state metrics using the 

branch metrics and the estimate of received symbols and; 
15 performing a backward recursion by calculating backward state metrics using 

the branch metrics and the estimate of received symbols stored from the forward recursion; and 
calculating a posteriori log-likelihood ratio of each select bit c in the estimate 

of received modulation symbols using the calculated forward state metric and the calculated 

backward state metric. 
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